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•/^ ! Abstract 

We investigate the internal space of Bessel functions which is associated to the 
group Z of positive and negative integers defining their orders. As a result we propose 
Qh! and prove a new unifying formula ( to be added to the huge literature on Bessel 

^ ' functions ) generating Bessel functions of real orders out of integer order one's. The 



unifying formula is expected to be of great use in applied mathematics. Some appli- 
cations of the formula are given for illustration 
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Very early studies [e.g., [0],0] proposed a unifying scheme for special functions sliown- 
ing that some of these functions may originate from the same structure .For Bessel func- 
tons of concern here their generating functions are representation " states" of the derivative 
and /or integral operators of arbitrary orders. More precisely we have the "inner" structure 
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Where / dz is the " truncated " primitive ,ie in defining the integral we omit the constant 
of integration / ^dz = f . For the polynomials such as Hermite and Laguerre for instance 
,the generating functions only involve the realization of the set N of positive integers, with 
slight modifications of the derivative operators to account for the conventions used in 
defining these polynomials.lt is important to note that although this common structure 
only set up the z-dependance of the generating functions ,it is the "dynamical" part of 
the scheme so to say. The t dependence is simply set by imposing some given desired 
properties. For Bessel functions we require a "sjnnmetry" between positive and negative 
indices that is J_„ = (— 1)"J„ ,while for the polynomials it is the natural property of 
orthonormality that is invoked. In this paper we unify Bessel functions of integer orders 
with those of real orders as a complement to the unification of special functions initiated in 
the above cited papers and subsequent papers. To our knowledge there is no such relation , 
between integer and real orders , known to exist. Both functions are usually independently 
defined either by their own differential equations or by their expansions in the form of a 
series or integral .The mechanism underlying the unification we propose is of a different 
nature than the one described above, but makes however use of the operators dm meZ 
which served to define Bessel functions of integer orders. 

Let us first introduce the method in a simple example to see how the mechanism works 
to convert an integer into a real. Suppose we are given an abstract state | n) neZ and a 
set of raising n(m), m > and lowering m < operators .Then it is easy to show , given 
that data that the state | n + A > is related to the state | n > through the following 
formula ( Provided the involved series converges which is the case here ). 

Il(m) 

|n + A>=exp(-A ^ (-1)'"^-^) | n > (2) 



meZ/(0) 



Fourier transforming the | n > state as 



n >= r ^e^"^ I ^ > (3) 



Where the 11 (m) operators acts on the | 6 > state by simple multiplication by the 
factor e*"^ .We get 
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In deriving the last line use have been made of the known formula X]m=i(~l)"^ ^^^^ — 

It is essential that the | n > states have the appropriate weights for the Fourier 
components otherwise the series defining the transformation simply diverges. Reduced 
Bessel functions (pniz) = fit into the scheme .We indeed have a set of raising and 
lowering operators 11 (m) = (— meZ |Q 

(-l)'"p;p0n(^) = Ct>n+m{z) UeZ (5) 

and the involved series does converge as it has almost the same structure as for the 
example above. We then propose and prove the unifying formula. 

= exp(-A E (-1)-^^) ^ (6) 

The opertors act on the as 

n(m)0„ = (t)n+m (7) 

To check the above formula ,we could naively express the Bessel functions 0n+m as 
entire series and then perform the sum over m .This is however a cumbersome procedure. A 
simple and illuminating way to proceed is to use the integral representation of 0„ , the 
"analog "of the | n > state Fourier decomposition. 
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Where r is a complex variable and Iq is a positively oriented closed path encircling the 
origin one time .Expanding the exponential in |^ as . 

Y\(m) 
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= t^-^i E i-ir^YMz). (9) 

p=0 P- meZ/iO) 

And using the propery that n(mi) Il{mp) = n(mi + mp), the term of order 

takes the form. 

;— 2^ mp (10) 



Inserting the 0„ integral representation into ITOl we get . 
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The relevant term to sum is . 

E (-ir^ (12) 

It is to be noted that the above series diverges on the whole complex plane except on the 
circle centring the origin and with half unity radius, on which it converges uniformally.lt 
is only on that circle do we have the right to commute the signs and / when we insert 
1^ into |TD| .Now the miracle happens .To make the series converge we just deform the path 
Iq down to the above circle as the integrand in |^ has only essential singularities at r = 
and r = oo. To compute the relevant sum ,put 2r = e'^ ,the series then converges to . 

meZ/{0) ^ ra=\ ^ 

= iO -n < 9 <n (13) 
= ln2T 

Where the branch cut of the logarithm is taken along the negative real axis as 6 ranges 



from —71 to TT .Inserting the result into |Tl| ,the latter becomes. 



Where the new positively oriented path 1 is now getting round the cut lying on the negative 



real axis. The above expression in 14 is however nothing than the order term of the 



real order Bessel function 0n+A when it is written in terms of its integral representation. 
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<Pn-,x{z) = — (-) r-("+^)- - -)rfr (15) 
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To see this ,rewrite (2r)^^ = exp(— A/n2r) and take the order which is ^(— A)^(/n2r)P 
.We thus have succeed to prove the proposed unifying formula. 

Now we propose two applications of this formula for illustration . The first immediate 
application is this .Bessel functions enjoy the essential property 

(-l)'"p;p0n(2:) = 0„+„(^) neZ (16) 
This formula is easily generalized to 0„+a with A real since the operator -^^^^ appearing 
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that we know the generating function of the reduced integer order Bessel function .That 
is we want to indirectly compute 

n=oo 

E <Pn+x{z)r (17) 
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Where t is a complex parameter .To perform the sum we use the unifying formula. 
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Where $(^,t) is the generating function for Bessel functions of integer orders. That is 
$(z,t) = Yl,'^=-QQ'Pn{z)t'^ .The action of the n(m) operators are simply . 

n(m)$ = t-™<l> (19) 



Then the third line in |I8| involves the following sum we got used to. 
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This series is only convergent on the circle of radius unity ,so that we put t = e*^ and 
find. 
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Putting this into [T^ we get the end result that the generating function for as defined 
above is shown to be t'^ times the generating function for Bessel functions of integer 
orders. Our formula shows that this is true on the cirle of radius unity except on the the 
point -1 . 
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E <Pn+x{z)e = t-^ E 0n(^)t" = t-'exp(---). (22) 
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Note that our formula only predict the result on the circle .If we assume in addition 
that the series is regular on the complex plane except on the branch cut on the negative 
axis, then the result we found is valid outside the circle and cover the whole domain of reg- 



a direct method .We will use the integral representation of the reduced Bessel function 
jOtherwise the computation is simply awful. 



n=oo J.-A n=oo „ j. 2 

n=-oo 27r«„_^7/ 2r 4r r 

Making the change of variable ^ = e*^.We rewrite the above expression as 

n=oo j-A "=oo -/- ~2 

E 0n+A(^)r = E e-^e^^^exp(^e-^^ - |^e^^))rf^ (24) 
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And knowing that the involved series converges to the sum of Dirac distributions. 

n=oo p=oo 

J2 e^"^= 27r(5(^ - 27rp) (25) 
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We get after insertion of the above result . 
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E 0n+A(^)t" = 6{e-27Tp)e^''exp{-e~''~-e^')de (26) 
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= t"^exp(--— ) 
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Note that only the p=0 term in the above sum contributes since the path of integration 
encircles the origin only one time. Here again we get the same result with the advantage 
however that the result is valid on the whole complex t plane except the cut .This is 
because we can always deform the r contour so to make I ^ |= 1 and hence to sum the 
series. Let us end up with the following worth remarks. First ,we have shown that Bessel 
functions of real orders can be obtained from integer orders ones through the mechanism 
we described above, and from this point of view integer order Bessel functions are more 
"fundamental" .Second,there is some very known functions which appear in mathematical 
physics but are simply constructed from various combinations of the linairely independant 
Jp{z) and J-p{z) , where p is real .These are Hankel Hp{z), Hp{z) and Neumann functions 
Np{z) which are expressed as 
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H\z) = Jp{z)+iNp{z) (27) 
H\z) = J,{z)-iNp{z) 

These functions are also concerned with the above unification but to give the formula 
which relates integer orders to real one's is not a straigthforward matter although they 
are simply linear combinations of Bessel functions. The reason is that the unifying formula 
acts on Jn{z)/ z^ and not directly on Jn{z) and therefore we cannot naively apply it to the 
defining formulas [2^ with p integer .We managed however to unify Neumann's functions 



of different orders using our unifying formula but following an indirect way .The result 
will appear in a forthcoming paper. 
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